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Abstract 

We extend the validity of Brill's axisymmetric positive energy the- 
orem to all asymptotically flat initial data sets with positive scalar 
curvature on simply connected manifolds. 
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1 Introduction 

In [3] Brill proved a positive energy theorem for a certain class of maximal, 
axi- symmetric initial data sets on M^. Brill's analysis has been extended 
independently by Moncrief (unpublished), Dain (unpublished), and Gibbons 
and Holzegel [9] to the following class of metrics: 

g = e-2C/+2a (^^2 ^ ^^2^ ^ p2^-2U ^ ^^^^^ ^ ^^^^)2 _ ^^ -^^ 

All the functions are assumed to be (/^-independent. 

The above form of the metric, together with Brill's formula for the mass, 
are the starting points of the recent work of Dain [7], who proves an upper 
bound for angular-momentum in terms of the mass for a class of maximal, 
vacuum, axi-symmetric initial data sets with a metric of the form above. 

The aim of this series of papers is to extend the validity of Brill's positiv- 
ity theorem, as well as that of Dain's inequality, to all maximal, asymptoti- 
cally flat, vacuum initial data sets (M, g, K) invariant under a U(l) action 
on simply connected manifolds. In fact, our analysis extends Brill's proof 
of positivity of mass to the above class of initial data sets, except that the 
condition of being vacuum is weakened to the requirement of positivity of 
the matter density. 

More precisely, in this paper we prove that any sufficiently differentiable, 
asymptotically flat, axially symmetric metric on can be written in the 
form (1.1). In general the functions appearing in (1.1) will not satisfy the 
fall-off conditions imposed in [7,9], but we verify that the proof extends to 
the more general situation. The result is further extended to include metrics 
with several asymptotically flat ends provided the manifold is simply con- 
nected. In the second paper of this series [6] the constructions of the current 
paper will be used to extend the validity of Dain's angular momentum in- 
equality to the class of metrics considered here. We will further allow those 
non- vacuum models which admit a twist potential to, see [6] for details. 

It is conceivable that, regardless of simple-connectedness and isotropy 
conditions, axi-symmetric metrics on manifolds obtained by blowing-up a 
finite number of points in a compact manifold can be represented as in 
(1.1), with the coordinates {p, z) ranging over a subset $7 of M^, and with 
identifications on dO,, but this remains to be seen; in any case it is not clear 
how to adapt the arguments leading to the mass and angular momentum 
inequalities to such situations. 
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2 Axi-symmetric metrics on simply connected asymp- 
totically flat three dimensional manifolds 

Let us start with a general discussion of Riemannian manifolds (M, g) with 
a Killing vector r/ with periodic orbits; without loss of generality we can 
assume that the period of principal orbits is 27r. 

Let M/U(l) denote the collection of the orbits of the group of isometries 
generated by r], and let vr : M — > M/U(l) be the canonical projection. An 
orbit p E M/U(l) will be called non- degenerate if it is not a point in M. 
Recall that near any p £ M/U(l) which lifts to an orbit of principal type 
there exists a canonical metric q defined as follows: let X,Y £ Tp(M/U(l)), 
let p £ M be any point such that irp = p, and let X ,Y £ TpM be the unique 
vectors orthogonal to r] such that 'it^:X = X and tt^:Y = Y. Then 

q{X,Y):=g{X,Y) . (2.1) 

(The reader will easily check that the right-hand-side of (2.1) is independent 
of the choice of p G iT^^{{p}).) 

There exists an open dense set of the quotient manifold A//U(l) which 
can, at least locally, be conveniently modeled on smooth submanifolds (per- 
haps with boundary), say A^, of M, which meet orbits of t] precisely once; 
these are called cross-sections of the group action. (For metrics of the form 
(1.1) there actually exists a global cross-section N, meeting all orbits pre- 
cisely once.) The manifold structure of Af/U(l) near p is then, by definition, 
the one arising from N. For 

p(z N := N\{r] = 0} 

and for X,Y e TpN set 

,(x.r) = »(x.n-^*ii4M^. (2.2) 

One easily checks that this coincides with our previous definition of q. 

The advantage of (2.2) is that it allows us to read-off properties of q 
directly from those of g near N. On the other hand, the abstract definition 
(2.1) makes clear the Riemannian character of q, and does not require any 
specific transverse submanifold. This allows to use different A^'s, adapted to 
different problems at hand, to draw conclusions about M /U(l); this freedom 
will be made use of in what follows. 
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Clearly all the information about g is contained in q and in the one-form 
field 

■=9{V,-) , 

since we can invert (2.2) using the formula, valid for any X,Y ^ TAI, 

,(x,y) = ,(p,x,p,y) + ^(!^4M^, (2.3) 

where P^j : TM TN is the projection from TM to TN along ij. (Recall 
that is defined as follows: since r] is transverse to TN, every vector 
X £ TM can be uniquely written as X = arj + Y, where Y £ TN, then one 
sets Pr^X := Y.) In order to establish (2.3) note, first, that this is only a 
rewriting of (2.2) when both X and Y are tangent to N. Next, (2.3) is an 
identity if either X or y is proportional to r], and the result easily follows. 

Let x^, A = 1,2 he any local coordinates on N, propagate them off N 
by requiring that ^r]X^ = 0, and let ip he & coordinate that vanishes on N 
and satisfies ^r^f = 1- Then rj = d^, and Prj{X^dA + X'^d^) = X^Oa, so 
that (2.3) can be rewritten as 

g = QAB dx^dx^ +g{r], rj) {dip + dAdx^ f , (2.4) 



with 

d^qAB = d^0A = d^{g{i],r])) = 



2.1 Global considerations 

So far our considerations were completely general, but local. Suppose, how- 
ever, that M is simply connected, with or without boundary, and satisfies the 
usual condition that it is the union of a compact set and of a finite number 
of asymptotically fiat ends. Then every asymptotic end can be compactified 
by adding a point, with the action of U{1) extending to the compactified 
manifold in the obvious way. Similarly every boundary component has to 
be a sphere [10, Lemma 4.9], which can be filled in by a ball, with the action 
of U{1) extending in the obvious way, reducing the analysis of the group 
action to the boundaryless case. Existence of asymptotically fiat regions 
implies (see, e.g., [2]) that the set of fixed points of the action is non-empty. 
It is then shown in [13] that, after the addition of a ball to every bound- 
ary component if necessary, M is homeomorphic to M'^, with the action of 
U{1) conjugate, by a homeomorphism, to the usual rotations of M^. On the 
other hand, it is shown in [12] that the actions are classified, up to smooth 
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conjugation, by topological invariants. It follows that the action is in fact 
smoothly conjugate to the usual rotations of M^. In particular there exists 
a global cross-section N for the action of U{1) away from the set of fixed 
points =e/, with N diffeomorphic to an open half-plane, with all isotropy 
groups trivial or equal to U{1), and with ^ diffeomorphic to M.^ 

Somewhat more generally, the above analysis applies whenever M can 
be compactified by adding a finite number of points or balls. A nontrivial 
example is provided by manifolds with a finite number of asymptotically flat 
and asymptotically cylindrical ends, as is the case for the Cauchy surfaces 
for the domain of outer communication of the extreme Kerr solution. 

2.2 Regularity at the axis 

In the coordinates of (1.1) the rotation axis 

£^ ■={g{r],7]) =0} 

corresponds to the set p = 0, which for asymptotically flat metrics is never 
empty, see, e.g., the proof of Proposition 2.4 in [2]. 

In order to study the properties of q near £/ /U(l) ~ £/, recall that £/ is 
a geodesic in M. It is convenient to introduce normal coordinates {x, y, z) : 

^ defined on an open neighborhood of a/, where z is a unit- 
normalized affine parameter on £/, and (x, y) are geodesic coordinates on 
exp((T^)-'-). Without loss of generality we can assume that ^ is invariant 
under the flow of rj. 

As is well known, we have (recalling that orbits of principal type form 
an open and dense set of M, as well as our normalization of 27r-periodicity 
of the principal orbits) 

rj = xdy- ydx . 

If we denote by (pt the flow of ry, on ^ the map (pT^ is therefore the symmetry 
across the axis ^ : 

(j)^{x,y,z) = {-X, -y,z) . 

This formula has several useful consequences. First, it follows that the 
manifold with boundary 

N := {x>0,y = 0} 

is a cross-section for the action of U(l) on ^ . This shows that near zeros 
of 7] the quotient space M/U(l) can be equipped with the structure of a 

am grateful to Joao Costa and Allen Hatcher for discussions or comments on the 
classification of U{1) actions. 
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smooth manifold with boundary. The analysis of the behavior of q near 
dN ~ £^ requires some work because of the factor 1/(7(7/, r/) appearing in 
(2.2). 

For further use we note that the manifold 

iV := {y = 0} C (2.5) 

provides, near =2/, a natural doubling of N across its boundary . 

In order to understand the smoothness of g on and N , we start by 
considering the function 

f{x,z) := c/(77,r?)(x,0,z) . 

Then f(—x,z) = f{x,z) because g{'i],rj) o 0^ = g{i],rj). It follows that all 
odd x-derivatives of / vanish at x = 0. It is then standard to show, using 
Borel's summation lemma (c/., e.g., [5, Proposition CI, Appendix C]), that 
there exists a smooth function h{s, z) such that 

f{x, z) = x^h{x'^, z) . 



Letting p = y^x^ + y^, invariance of g under (f)t allows us to conclude that 
g{ri, r]){x,y, z) = g{r], r]){p, 0, z) = p^h{p^ , z) . (2.6) 
Define (p via the equations 

X = p cos if , y = p sin ifi , 

so that 

V = d0 . 

Considerations similar to those leading to (2.6) (see Lemma 5.1 of [5]) show 
that there exist functions a, (3, 7, 5, p, and gzz, which are smooth with respect 
to the arguments and z^ with 

;u(0,z) = l, c/i5(0,i) = l, 

such that 



g = gad-z + 2apdzdp + 2l3p dzdip + ^p dp + 25p dpd(p + p{dp + p d(p 

=2 ]az +2[a ] pazap + ( n + jjj 

p. / \ p / \ p 



6 /3 \2 

+up^[d(2i + -pdp+ -dz . (2.7) 
' p p J 



^By this we mean that a(s, z) is a smooth function of its arguments, and enters (2.7) 
in the form a(p^, z), etc. 
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We say that is a doubling of a manifold N across a boundary if 
consists of two copies of N with points on identified in the obvious way. 
From what has been said, by inspection of (2.7) it follows that: 

Proposition 2.1 The quotient space M/U(l) has a natural structure of 
manifold with boundary near si . The metric q and the one- form 9 are 
smooth up-to-boundary, and extend smoothly across £/ by continuity to 
themselves when M/U(l) is doubled at £/ . 

For further use we note the formula 

giv,v)=P^ + 0{p^) , (2.8) 

for small p, which follows from (2.7), where p is either the geodesic distance 
from s^, or the geodesic distance from £/ on exp((T£/)-'-) (the latter being, 
for small p, the restriction to exp((T=e/)-'-) of the former). 

2.3 Asymptotic flatness 

We will consider Riemannian manifolds (Af , g) that are asymptotically flat, 
in the usual sense that there exists a region Mcxt C M diffeomorphic to 
\ B{R), where B{R) is a coordinate ball of radius R, such that in local 
coordinates on Mext obtained from \ S (i?) the metric satisfies the fall-off 
conditions, for some A; > 1, 

9ij - S^j = Ofc(r-V2) ^ (2.9) 
dkgij G L2(Mext) , (2.10) 
R'jki = o{r~^'^) , (2.11) 

where we write / = Ofc(r°) if / satisfies 

dk,...dkj = o{r''-') , 0<£<k. 

It is well known that (2.9)-(2.10) together with R{g) > or R{g) £ L^, 
where R{g) is the Ricci scalar of g, guarantees a well-defined ADM mass 
(perhaps infinite). On the other hand, the condition (2.11) (which follows 
in any case from (2.9) for /c > 2) is useful when analyzing the asymptotic 
behavior of Killing vector fields. 

We will use (2.9)-(2.11) to construct the coordinate system of (2.3), and 
also to derive the asymptotic behavior of the fields appearing in (2.3). We 
start by noting that the arguments of [1, Appendix C] with A^ = there 
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show that there exists a rotation matrix oo such that in local coordinates on 
Mext we have 

7?^ = u;>^>Ofc(ri/2), (2.12) 

where is anti-symmetric. It will be clear from the proof below (see 
(2.23)) that this equation provides the information needed in the region 

x2+y2>22, x2 + y2+22>^2_ ^2.13) 

However, near the axis a more precise result is required, and we continue by 
constructing new asymptotically flat coordinates which are better adapted 
to the problem at hand. The difficulties arise from the need to obtain decay 
estimates on g — 5, where 5 is the Euclidean metric on M^, and on 9, which 
are uniform in r up to the axis ^ . 

Let (i*) = {x,y,z) be coordinates on \ B{R), obtained by a rigid 
rotation of x*, such that uj^jx^ = yd^ — xdy. Set 

— xo_^ ^ y:=^ — y-^ , ^ •= ^ I o </.^ . (2.14) 

Using the techniques in [1, 2] one finds 

0s(x*) = (cos(s)2;— sin(s)y+z^(s, X*), sin(s)x+cos(s)y+2;^(s, X*), i+z^(s, X*)) , 
with satisfying 

We then have 

dz _ ^ 1 [^^ dz^{(ks{x')) 



OZ 2lT Jq oz 



Further, 



similarly 



dx 2tt Jq dx 



dy 

The estimates for the derivatives of x and y are straightforward, and we 
conclude that 
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where, by an abuse of notation, we write again for the functions (x, y, z). 
Standard considerations based on the imphcit function theorem show that, 
increasing R if necessary, the j;*'s form a coordinate system on \ B{R) in 
which (2.9)-(2.11) hold. Subsequently, (2.12) holds again. 
From (2.14) one clearly has 

VsGM z o = z , 

which shows that the planes 

^^■={z = t} , t€R, \t\>R, 

are invariant under the flow of 77; equivalently, 

r/^ = . 

Moreover, 

XO(j)^ = -X, yo(j)^ = -y, (2.15) 

SO that all points with coordinates x = y = are fixed points of 0,^, and 
that these are the only such points in M^xt- Equation (2.15) further implies 
that (pT^ maps the surfaces {x = 0} and {y = 0} into themselves. Since 0,^ 
is an isometry, we obtain 

9ab{0, y, z) = gab{0, -y, z) , g^ziO, y, z) = gzz{0, -y, z) , 

9za{0, y, z) = -gza{0, -y, z) ; (2.16) 

similarly 

gab{x, 0, z) = gab{-x, 0, z) , gzzix, 0, z) = gzz{-x, 0, z) , 

gza{x, 0, z) = -gzai-x, 0, z) . (2.17) 

Equation (2.16) leads to 
^^(0,0,.)=0, ^^(0,0,z)=0, ^(0,0,z)=0 (2.18) 

for ^ G N (or at least as far as the differentiability of the metric allows) . The 
analogous implication of (2.17) allows us to conclude that 

^(0,0,z)=0, ^(0,0,z)=0, gaz{0,0,z)=0. (2.19) 
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Incidentally, the last two equations in (2.19) show that {x = y = 0} is a 
geodesic; this follows in any case from the well-known fact that the set of 
fixed points of an isometry is totally geodesic. 

Consider a point p lying on the axis of rotation £/, then (j)tip) = P ^or 
all t, in particular (pT^ip) = P- From what has been said we obtain that 

£/ n Mext C{x = y = 0} . (2.20) 

Recall, again, that every connected component of the axis of rotation £/ 
is an inextendible geodesic in {M,g). Since the set at the right-hand-side 
of (2.20) is a geodesic segment, we conclude that equality holds in (2.20). 
Hence 

7]\0,0,z)=0 (2.21) 

and, for \z\ > R, the origin is the only point within the plane at which 
1] vanishes. 

We are ready now to pass to the problem at hand, namely an asymptotic 
analysis of the fields g{r],rj), q and 9 as in (2.4); we start with q. For p 
sufficiently large the hypersurface {y = 0} is transverse to rj (for small p we 
will return to this issue shortly) and therefore the coordinates 

(x^) := ix,z) 

on this hypersurface, with x > 0, can be used as local coordinates on 
M/U(l). The contribution of gAB to qAB is of the form qab = ^AB + 
Ofc(''~^^^)) which is manifestly asymptotically flat in the usual sense. Next, 
from (2.9) and (2.12) we obtain 

5(?7,?7) =p2 + Ofc(r3/2) ; (2.22) 
here, as elsewhere, p'^ = x'^ -\- y'^. Further 



gMV 9B,V ^^A^^B ^ {5M + ou{r-'l')){u:\x'^ + o,{r''^)) 

5Bj + Ok{r-^l^))Uhx'' + ou{r^l^) 



P^+Ok{r^/^) 
Ok{r^^'^)dx^dx^ 
p"^ + Ok{r^/'^) 



dx^dx^ 
(2.23) 



because u!^aX°'^Hx''dx^dx^ = {xdy — ydx)'^ , which vanishes when pulled-back 
to {y = 0}. In the region (2.13) we thus obtain 

QAB = Sab + Ok{r^^/^) , (2.24) 
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which is the desired estimate. However, near the zeros of rj this calculation 
is not precise enough to obtain uniform estimates on q and its derivatives. 

In fact, it will be seen in the remainder of the proof that we need uniform 
estimates for derivatives up to second order. Since g{r], r]) vanishes quadrat- 
ically at the origin we need uniform control of the numerator of (2.23) up to 
terms 0{p^), in a form which allows the division to be performed without 
losing uniformity. 

So in the region {p < \z\} n M^xt, in which \z\ is comparable with r, we 
proceed as follows: Let 

X\ = X\{z) := 1^(0, 0, z) , Xab := 9ac{0, 0, z)A'=fe ; 

note that A"^;, = lo"-}, + Ofc_i(|z|~^/^) = u;"-h + Ofc_i(r~^/^), similarly for Xab- 
The Killing equations imply that A^;, is anti-symmetric, hence 

Kx = Xyy = , Xxy = -Xyx = 1 + Ofc- 1 ( | Z | " "^^^ ) = 1 + Ofc_l(r~^/^) . 

From (2.21) we further obtain 

Recall the well known consequence of the Killing equations, 
which implies, at £/, 

= VaVbVc = dadbVc , (2.25) 

= VaVbVz = daVbT], - T'azXbc = OadbVz " 2r^„,Afee • (2.26) 

From (2.12) we obtain, by integration of third derivatives of i]a along rays 
from the origin x = y = within the planes z = const, 

^'"^'^ =o._3(k|-^/^)x^ = o,_3(r-^/^)x% 



dx^dx' 

and then successive such integrations give 

Va = Xabx' + Ofc_3(r-5/2)^c^'^^e _ ^2.27) 
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At 2/ = we conclude that 

Similarly we have jj]^ = R^ij^rjg, hence Va^tv'^ = dadh't]'^ = at , 
and we conclude that 

r]" = X\x^ + Ofc_3(r-5/2)^c^d^e _ (2.28) 

This allows us to prove transversality of rj to the plane {y = 0}. Indeed, 
from (2.28) at y = we have 

7f = (1 + o{r~^/^))x + o(r-5/2)^3 ^ ^ o{r-^/^))x 

which has no zeros for x ^ and r > if i? is large enough. Recall that we 
have been assuming that |x| < \z\ in the current calculation; however, we 
already know that rj is transverse for \z\ > and transversality follows. 
Increasing the value of the radius R defining Mext if necessary, we conclude 
that {y = 0, X > 0} n Mext provides a global cross-section for the action of 
U(l) in Mext. 

Using (2.26), a similar analysis of rjz gives 

Vz = - r^azU Ab,xV + Ofc„3(r-5/2)^c^d^e ^ 
Ofc-i(r-3/2) 

We are now ready to return to (2.22), 

giv, V) = mv' = VaV" = + Ofc-3(r^^/')x"x''x^x'^ , (2.29) 
where, at y = 0, 

:= gatWx'X'dx'' = (1 + Ofc_i(r-V2))^2 . 

it follows that the last equality also holds for g{r],ri) with k — 1 replaced by 
k — 3. Instead of (2.23) we write 

gAiV'gBjV^ ^^A^^B _ r]Ar]Bdx^dx^ 



g{r],rj) (1 + Ofc_3(r-V2))a;2 

r/2 dx^ + 2rjxrizdx dz + rjl dz 



(l+o,„3(r-V2))x2 

Ofc_3(r~^)x2(ix^(ix^ 
(l+Ofc_3(r-i/2)) 

Ofc_3(r~^)dx^dx^ . (2.30) 
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We conclude that (2.24) holds throughout {y = 0} n Afcxt with k replaced 
by A; - 3. 

To analyse the fall-off of Bp and Az , note first that the discussion in the 
paragraph before (2.4) shows that it suffices to do this at one single surface 
transverse to the flow of the Killing vector field rj; unsurprisingly, we choose 

N ■={y = 0, x>0, + z2 > ^ 

with R sufficiently large to guarantee transversality. Next, from (1.1) we 
find 

r]idx' =g{r],-) =g{d^,-) = g{r],r]){dip + pBpdp + A^dz) , 

which will allow us to relate Bp and A^ to r]i if we determine, say dup and dip 
on N. For the sake of clarity of intermediate calculations it is convenient to 
denote by z the coordinate z appearing in (1.1), we thus seek a coordinate 
transformation 

{x, y, z) (p, if, z) , with z = z everywhere and p = x on N, 

which brings the metric to the form (1.1), with z there replaced by z. We 
wish to show that, on N, 








7]^ 




H 





if 














1 1 = . (2.31) 



The second column is immediate from 

dx_Q _^ dy_Q ^ dz_ 

dip ^ dip ^ dip 



ri^d, + rfdy + ij'd, =r] = d^ = —d, + ^dy + —d, . 



Similarly the third row follows immediately from dz = dz. It seems that 
the remaining entries require considering J^^ . Now, 93 is a coordinate that 
vanishes on A^, so that dxip = dz^p = there. From rfdip = 1 we thus obtain 
dyip = ^-/rf . Next, p = x on N, giving dxP = 1 and dzP = there. The 
equation ij'^dip = gives then r/^ + rfdyp = 0, so that dyp = —ij^/rj^. The 
derivatives of z are straightforward, leading to 



dx dy dz 

Inverting leads to (2.31) 



[ m ^ m \ /I -v^'/riy 

|£ ^ ^ = l/r]y 

ax ay oz \ \ ' ' 

91 as dz \ 1 
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Prom now on we drop the bar on z. From (2.31) one immediately has 

on 



Similarly, again on A^, 

^ PaiViV) dp xg{r],r]) 

Finally, we note that 
9{V,V) 



otherwise, 



(2.32) 



Ofc_3(r 5/2), |x| < \z\, 
Okif'^^"^), otherwise. 



(2.33) 



-2U 



P' 



l + Ofe_i(r i/2)+Ofc_3(r 5/2)x2, |x| < |z|, 
1 + Ofe(r^^/2), otherwise. 



= l+o,_3(r-i/2^ 
In summary: 



(2.34) 



Proposition 2.2 Under (2.9) with k > 3 the metric q is asymptotically 
flat. In fact, there exist coordinates {x,y,z) satisfying (2.9) and a constant 
R > such that the plane {y = 0} O {r > R} is transverse to rj except at 
X = z = where rj vanishes and, setting x^ = {x, z) we have 



Furthermore (2.32)- (2.34) hold. 



(2.35) 



2.4 Isothermal coordinates 

We will use the same symbol q for metric on the manifold obtained by 
doubling M/U(l) across the axis. 
We start by noting the following: 

Proposition 2.3 Let q be an asymptotically flat metric on M? in the sense 
of (2.35) with k > 5. Then q has a global representation 

q = e2«(di;2 + dw^) , with u . (2.36) 

In fact, u = Ofc_4(r~-'^/2). 
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Remark 2.4 The classical justification of the existence of global isothermal 
coordinates proceeds by constructing the coordinate v of (2.36) as a solution 
of the equation AgU = 0. A more careful version of the arguments in the 
spirit of [16, Lemma 2.3] shows that v has no critical points. However, the 
approach here appears to be simpler. 

Proof: Let qab = &~'^^(1aBi then q is flat if and only if u satisfies the 
equation 

A,n = , (2.37) 

where R{q) is the scalar curvature of q. For asymptotically flat metrics q, 
with asymptotically Euclidean coordinates [x^z), this equation always has 
a solution such that 

u+;uln(\/x2 + z2) — ^^jTT^^oo ' ^^e^e = 7- / Riq) dfiq , (2.38) 

where dfiq is the volume form of q. More precisely, we have the following: 

Lemma 2.5 Consider a metric q on satisfying 

QAB - ^AB = Of (r"^/^) 

for some £ > 2, with {x^) = {x,z). For any continuous function R = 
Of_2(r~^/^) there exists u = Of_i(r~^/^) and a solution of (2.37) such that 



u — ^ln(\/ x"^ + z^) 



with fi as in (2.38). 



Proof: We start by showing that (2.37) can be solved for |x| large. Indeed, 
consider the sequence Vi of solutions of (2.37) on the annulus 

T{p,p + i) :=D{0,p + i)\D{0,p) , 

with zero boundary values. Here p is a constant chosen large enough so 
that the functions ztClxl"-*^/^, with C = 8||i?|x|^/^||i=o, are sub- and super- 
solutions of (2.37). The usual elliptic estimates show that a subsequence 
can be chosen which converges, uniformly on compact sets, to a solution 
V = Of„i(r^V2) of (2.37) on \ D{0, p). In the notation of [4] we have in 
fact V G C^Sii2fi ^^^y ^ ^ Furthermore, using the techniques in [4] 

one checks that v = oe-i{r~^^'^). 
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We extend v in any way to a C^"^'"^ function on M^, still denoted by 
V. Let q := e~'^'"q, then q is flat for |x| > p. Let be any g-parallel 
orthonormal co-frame on \ D{0,p), performing a rigid rotation of the 
coordinates if necessary we will have = dx^ + 0£_i{\x\~^^'^)dx^ for 
I a:; I large. Let x^ be any solutions of the set of equations dx^ = e^. By the 
implicit function theorem the functions x"^ cover \ D{0,p), for some p, 
and form a coordinate system there, in which (jab = ^AB- 

Since (2.37) is conformally covariant, we have reduced the problem to 
one where R has compact support, and g is a C^~^''^ metric which is flat 
outside of a compact set. This will be assumed in what follows. 

Let us use the stereographic projection, say ip, to map to a sphere, 
then (2.37) becomes an equation for u := (u — v) o on \ {i^}, where 
is the north pole of S'^, of the form 

A,,n = |x|V, (2.39) 

where Hab ■= \x\-^qAB is a C'~^'^ metric on S\ similarly / is a C'~^ 
function on S'^ supported away from the north pole. In fact, in a coordinate 
system 

/ = (2.40) 

near = {y^ = 0}, where the x^^s are the explicitly flat coordinates on 
\ D{0, R) for the metric q, we have 

hAB = Sab ■ 

Let i^fc(S'^) be the usual L^-type Sobolev space of functions on S'^ and 

set 

Hfc = {x G Hk{S^) I I^^Xdfih = O} , (2.41) 
where dp,h is the measure associated with the metric h. We have 

Proposition 2.6 Let h be a twice- differentiable metric on S"^, then A/j : 
M2 — > Ho is an isomorphism. 

Proof: Injectivity is straightforward. To show surjectivity, let X C 
be the image of BI2 by Ah, by elliptic estimates X is a closed subspace of 
L^(S'^). Let 99 G be orthogonal to X, then 

VXGM2 I^AhXdfih = 0. 
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Thus 99 is a weak solution of A/,,(/5 = 0, by elliptic estimates ip £ M2- But 
setting X = V integrating by parts one obtains dip = 0, hence (p is 
constant, which shows that X = Hq. □ 

Returning to the proof of Lemma 2.5, we have seen that (2.37) can be 
reduced to solving the problem 

A-^u = f, (2.42) 
where h is flat outside of a compact set. Let 

then 



Jr2 ^ ^ P-^"^ Jc{0,p) ^ 

= 27rfi = - fdfi-^ 



n 



Thus (2.42) is equivalent to the following equation for the function u :- 



and the right-hand-side has vanishing average. Transforming to a problem 
on 5^ as in (2.39), we can solve the resulting equation by Proposition 2.6. 
Transforming back to M? , and shifting u by a constant if necessary, the result 
follows. □ 

Returning to the proof of Proposition 2.3, we claim that /U = 0; that is, 

/ Riq)dfig = 0. (2.43) 

This is the simplest version of the Gauss-Bonnet theorem, we give the proof 
for completeness: consider any metric on satisfying 

QAB - Sab = 01(1) , R{q) G . 

Let O"", a = 1,2, be an orthonormal co-frame for q obtained by a Gram- 
Schmidt procedure starting from (dx^, dx"^), with connection coefficients w'^^. 
Then cj'^t = o(r^^). It is well known that, in dimension two, 

R{q) dfig = 2duj^2 ■ (2.44) 
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Equation (2.43) immediately follows by integration on B{R), using Stokes' 
theorem, and passing to the limit i? — > oo. 

To finish the proof, note that the metric qis a complete flat metric on , 
and the Hadamard-Cartan theorem shows the existence of global manifestly 
flat coordinates, say (vjw) so that q can be written as in (2.36). □ 

Returning to the problem at hand, recall that the metric q on has 
been obtained by doubling M/U(l) across £/. Let us denote by (j) the 
corresponding isometry; note that in Mext/U(l), in the coordinates {x,z) 
constructed in Section 2.3, the isometry cj) is the symmetry around the z- 
axis: </>(x,2;) = {—x,z). Similarly, in geodesic coordinates centred on £/, 
4>{x,z) = {-x,z). 

As (p is an isometry of q, preserving the boundary conditions satisfied by 
u, uniqueness of solutions of (2.37) implies that uo(j) = u. Smoothness on the 
doubled manifold shows that on £/ the gradient Vn has only components 
tangential to £/. This implies that is totally geodesic both for q and q. 

Choose any point p on £/. By a shift of {v, w) we can arrange to have 
{v{p),w{p)) = (0,0). Let {p,z) be coordinates obtained by a rigid rotation 
of {v, w) around the origin so that the vector tangent to =2/ at p coincides 
with dz- Then the axis {{0, z)}z£M. is a geodesic of q, sharing a common 
direction at p with £/, hence 

£/ = {{0,z)}zeR ■ 

Since (p is an isometry of q which is the identity on a/, it easily follows that 

so that M/U(l) = {p > 0}. We have thus obtained the representation (1.1) 
of g. 

The reader might have noticed that the function u constructed in this 
section is a solution of a Neumann problem with vanishing data on the axis. 

For further use, we note that from (1.1), on exp((Ti2/)-'-) the geodesic 
distance p from the origin equals 

and comparing with (2.8) we obtain 

q(0,z)=0. (2.45) 

Now, the function u = Ofc_4(r~^/^) of Proposition 2.3 equals u = 2(a — 
U) (compare (1.1)). By (2.45) and an analysis of Taylor expansions as in 
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Section 2.3 we infer that, at {y = 0}, 

a = Ofc_5(r-^/2^x . (2.46) 
From Proposition 2.2 we conclude: 

Theorem 2.7 Let k > 5. Any Riemannian metric on M.^ invariant under 
rotations around a coordinate axis and satisfying 

g^, - S,, = Ok{r-^/^) (2.47) 

admits a global representation of the form (1-1), with the functions U, a, 
Bp and satisfying 

^, = Ofc_3(r-3/2); Bp = Ok.s{r-'/^); U = Oksir-'/^) ; a = Ok-4{r-'/') . 

(2.48) 

Furthermore (2.46) holds. 

Remark 2.8 The decay rate o(r~^/^) in (2.47) has been tailored to the 
requirement of a well-defined ADM mass; the result remains true with decay 
rates o(r~") or 0{r~°') for any a G (0, 1), with the decay rate carrying over 
to the functions appearing in (1.1) in the obvious way, as in (2.48). 

2.4.1 Several asymptotically flat ends 

The above considerations generalize to several asymptotically flat ends: 

Theorem 2.9 Let k > 5, and consider a simply connected three-dimensional 
Riemannian manifold (M, g) which is the union of a compact set and of N 
asymptotically flat ends, and let Mgxt denote the first such end. If g is 
invariant under an action o/U(l), then g admits a global representation of 
the form (1.1), where the coordinates {z,p) cover (M x M+) \ {aj}^2; ^^^^ 
the punctures Si = (0, a^) lying on the z-axis, each Si representing "a point 
at infinity" of the remaining asymptotically flat regions. The functions U, 
a. Bp and A^ satisfy (2.48) in Mgxt- 
// we set 

then we have the following asymptotic behavior near each of the punctures 

U = 2\nri + 0k-A{r]'^) , « = Ok-M'^) , (2.49) 

where f = oi{r\^'^) means that ■ ■ -dAjf = Oi-j{ry^ ■') for < j < £. 
Finally, (2.45) holds. 
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Proof: As discussed in Section 2.1, M is diffeomorphic to minus a finite 
set of points and, after perfoming a diffeomorphism if necessary, the action 
of the group is that by rotations around a coordinate axis of . As in the 
proof of Theorem 2.7 there exists a function v = Ofc_4(r~^/^) so that the 
metric e~'^'"q is flat for large enough in each of the asymptotic regions. 
Equation (2.37) is then equivalent to the following equation for u — v, 

A,-..,(n -v) = -e2- + A,v) , (2.50) 

where the right-hand-side is compactly supported on M /U(l). Let Mcxt/U(l) 
be the orbit space associated to the first asymptotically flat region and let 
^l' be any smooth strictly positive function on Af/U(l) which coincides with 
in each of the remaining asymptotically flat regions of Af/U(l), where 
the y^^s are the manifestly flat coordinates there, with i/^ equal to one in 
Mcxt/U(l). Then (2.50) is equivalent to 

A^,-2.,(n -v) = -V^-ie2- + A,^;) . (2.51) 

Both the metric tpe~'^'"q and the source term extend smoothly through the 
origins, say i^, j = 1, . . . , N , of each of the local coordinate systems := 
y^/\y\'^- Simple connectedness of the two-dimensional manifold 

JV := M/U(l) U 

implies that JV w R'^, so that (2.51) is an equation to which Lemma 2.5 
applies. We thus obtain a solution, say w, of (2.51), and subsequently a 
solution V + w of (2.37) which tends to a constant in each of the asymptot- 
ically flat regions (possibly different constants in different ends), except (as 
will be seen shortly) in Mext where it diverges logarithmically. Note that at 
large distances in each of the asymptotically flat regions the function w is 
harmonic with respect to the Euclidean metric, hence approaches its asymp- 
totic value as with gradient falling-off one order faster. Similarly v 
has controlled asymptotics there, as in the proof of Lemma 2.5. Integrating 
(2.37) over M/U(l) one finds that the coefficient of the logarithmic term is 
again as in (2.38). 

In order to determine that coefficient, we note that since KiR"^ there 
exists a global orthonormal coframe for g, e.g. obtained by a Gram-Schmidt 
procedure from a global trivialization of T*R^. As a starting point for this 
procedure one can, and we will do so, use a holonomic basis dx^ with the co- 
ordinate functions equal to the manifestly flat coordinates in Moxt/U(l). 
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Furthermore, after a rigid rotation of the y^^s if necessary, where the y^^s 
are the manifestly flat coordinates for the metric e~'^^'^^'"\ in the asymptot- 
ically flat regions other than Mext/U(l), we can also assume that the dx^'s 
coincide with d{y^/\y\'^) near each i'j. By (2.44) and by what is said in the 
paragraph following that equation we have 

N 

where the C(i^,e)'s are circles of radius e centred at the z^'s. Near each 
the metric q takes the form e^^'"+'^^ 5ABdy^dy^ = e^^'"+'^^x\~^6ABdx^dx^ . 
The co-frame 9^ is given by 6^ = e^'"'^'^^ \x\~'^ dx^ , leading to 



Co\ = -^{X^dx^ - X^dx^) + 0(|f|-^/2)^^A 

\x\ 



so that 



lim (h uj^2 = 'i'Tr . 



We note that we have proved: 

Proposition 2.10 Let q be a Riemannian metric on a simply connected 
two-dimensional manifold which is the union of a compact set and N ends 
which are asymptotically fiat in the sense of (2.35), then 



4. 



]^ J R{q)dfi, = 2{N - 1) 



□ 



Since /U / 0, the function v + w obtained so far needs to be modified to 
get rid of the logarithmic divergence. In order to do that for j = 2, . . . ,N 
we construct functions Uj, q-harmonic on M/U(l), such that, in coordinates 
x^ which are manifestly conformally flat in each of the asymptotic regions, 

r ln|x|+o(l), inMe,t/U(l); 
Uj = < — In -|- 0(1), in the J 'th asymptotic region; (2.52) 
0(1), in the remaining asymptotic regions. 

This can be done as follows: let uj be any smooth function which in lo- 
cal manifestly conformally flat coordinates both near and on Mext/U(l) 
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equals In | x | , and which equals one at large distances in the remaining asymp- 
totically flat regions. Let ip be as in (2.51), then A^^-2(v+w)qUj is compactly 
supported in M/U(l). Further 



JM/Vd) 



M/U(l) 

lim / In |x| • n — lim / Dlnlxj-n 

. 



We can therefore invoke Lemma 2.5 to conclude that there exists a uni- 
formly bounded function v, approaching zero as one recedes to infinity in 
Mext/U(l), such that 

A^^-2iv+vj)gV = -A^^-2{v+-w)gUj . 

Subsequently the function uj := uj + v is g-harmonic and satisfies (2.52). 
The function 

N 

u := V + w + 2 Uj + a , 

where a is an appropriately chosen constant (compare [4]), defines the de- 
sired conformal factor approaching one as one tends to infinity in Mcxt/U(l) 
so that e~^"(7 is flat. This conformal factor further compactifies each of the 
asymptotic infinities except the first one to a point, so that e~'^^q extends 
by continuity to a flat complete metric on the simply connected manifold 
ryV. By the Hadamard-Cartan theorem there exists on ^ a global manifestly 
flat coordinate system for e"^"(7. The axis of rotation can be made to co- 
incide with a coordinate axis as in the proof of Theorem 2.7. It should be 
clear that the points at infinity lie on that axis. 

In order to prove (2.49), note that the construction above gives directly. 

1 /2 

U - a = u = Ci + 2lnri + Ofc„4(r • ) , 
Next, U can be determined by applying an inversion 

y^^{p,z-ai) = {x^) = iy^/\y\^) (2.53) 

to (2.34), 
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Since a vanishes on the axis (y^)^ + (y^)^ = in each of the asymptotic 
regions, we conclude that Cj = 0, and (2.49) fohows. □ 



3 ADM mass 

Let m be the ADM mass of g, 



m :- 



where dSi = di\{dx A dy A dz). This has to be calculated in a coordinate 
system satisfying (2.9). Typically one takes Sr to be a coordinate sphere 
S{R) of radius R; however, as is well-known, under (2.9) Sr can be taken 
to be any piecewise differentiable surface homologous to S{R) such that 

mi{r{p)\p £ Sr] ^r.^oo oo . (3.1) 

We will exploit this freedom in our calculation to follow. 

We introduce new coordinates x and y so that p and tp in (1.1) become 
the usual polar coordinates on R^: 

X = p cos if , y = p sin ip . 

This implies 

pdp = ^d{p'^) = xdx + ydy , 

p^dip = xdy — ydx , 
p'^dip^ = dx^ + dy"^ - dp^ . 

Inserting the above in (1.1) one obtains 

g = e-^''{dx^ + dy^)+^—^^^^^^{xdx + ydyf+e-^^+^''dz^ 

^ V ' P ^ V ' 

+2e~'^^{xdy - ydx) (^Bp{xdx + ydy) + A^dz^ 
+ terms quadratic in {Bp^Az) ■ (3-2) 
This will satisfy (2.9) if we assume that 

(e^" - l)x^ (e^° - l)xy (e^" - l)y^ _ 

Bpx'^ , Bpxy , Bpy^ , A^x , A^y = oxir^^l'^) , (3.4) 
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consistently with Theorem 2.7. Then the term occurring in the last line of 
(3.2) will not give any contribution to the mass integral: 



1 



g = e-"^^ {dx"^ + dy"^) + ^ — {xdx + ydyf + e 



(a) 



+ 2{xdy — ydx) (^Bp{xdx + ydy) + A^dz^ 



id) 



+01 (r ^)dx'dx^ . 



(3.5) 



Let us denote by the variables x,y. As an example, consider the contri- 
bution of (c) to the mass integrand: 



dSa 



(c) > gzz,zdSz - gzz,idSi = -gzz,adSa = (2{U - Q),a + o( 

A similar calculation of (a) easily leads to 

(a) + (c) {4Ui + o{r-^))dSi-2a,adSa . 



The contribution of (6) to the mass integrand looks rather uninviting at first 
sight: 



-,2a 



1 



+ 



xy 



-,2a 



1 



xy 



p2 



-,2a 



1 



X^ + 



dSr 



dSn, 



{x^ + y'')dSz . 



Fortunately, things simplify nicely if Sr is chosen to be the boundary of the 
solid cylinder 

Cr := {-R < z<R, 0< p< R} . (3.6) 

Then Sr is the union of the bottom Br = {z = —R , 0<p< R}, the lid 
Lr = {z = R , < p < R}, and the wall Wr = {-R < z < R , p = R}. 
On the bottom and on the lid we only have a contribution from dSz, which 
equals 

- (2a^z + o(^~^)) dx A dy 
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on the lid, and minus this expression on the bottom. On the wall dSz gives 
no contribution, while 

dSxlwn = {dy'^dz)\wii = x\wiidipAdz , dSylwn = —{dxAdz)\wfi = y\wiidipf\dz . 

Surprisingly, the terms in containing derivatives of a drop out, lead- 

ing to 



\Wr 



2a + o{r '^)^dipAdz 



We continue with the contribution of Bp to (d): 



((a;^ - y^)Bi^,y- {2xyBp,^x 

(2) 



dSx + 



dSy . 



(1) 



(3) 



(4) 



It only contributes on the wall Wr, giving however a zero contribution there: 
(x^ - y'^){xdy + ydx) + 2xy{ydy - xd^)) Bp dip A dz 

dip A dz = . 



(l)+{3) (4)+(2) 

{x^ + y^) {xdy - ydx)Bp 



Finally, produces the following boundary integrand: 
-ydzA;,dSx + xdzAzdSy + {xdy - ydx)Az 



dSz , 



and one easily checks that the dSx and dSy terms cancel out when integrated 
upon Wr, while giving no contribution on the bottom and the lid. 
Collecting all this we obtain 



m = lim 



: [ diUdSi + 2 [ 



X" 



4/ diUdSi + 2 I {a daa)d<pdz 

Wr P 



2 / dzOidxdy + 2 I dzCtdxdy 



lim — 



Br 

1 If 

di{U a)dSi-\ — / adipdz 

Sr '2 '2 Jwr 



We have the following formula for the Ricci scalar ^^^R of the metric (1.1) 
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(the details of the calculation can be found in [9]):^ 



„-2U+2a -1-1 1 r)rv, n2^>-2a 



(3.7) 

The Laplacian and the gradient D are taken with respect to the flat 
metric 6 on M^. 
Now, 



lim — 

R-^oo Air 



di{U-la)dSi + l 



a dip dz 



Wr 



lim 



1 /" r . ot. 1 da, o 1 , 

-y^JA,(C/--) + --ld^x+-/ .(p = 0..)<i. 



R 
R 



(3.8) 



The last integral vanishes by (2.45). Equations (3.7)-(3.8) and the domi- 
nated convergence theorem yield now 



m 



1 

16^ 



+— / (DUfd^x. 
Svr J 



(3.9) 



Since ^^^R = IGvr/i + KabK""^ > for initial data sets satisfying tr^K = 0, 
this proves positivity of mass for initial data sets as considered above. 
Suppose that m = with ^^^R > 0, then (3.9) gives 

= pBp^, - A,,p = DU = 0. (3.10) 

The last equality implies U = 0, and from (3.7) we conclude that 

1 da 

The maximum principle applied on the set 

B{R) \ {p < 1/R} 



'^In the time-symmetric case (3.7) can be viewed as a PDE for U given the remaining 
functions and the matter density. Assuming that this equation can indeed be solved, this 
allows us to prescribe freely the functions a, Bp and A^. In such a rough analysis there 
does not seem to be any constraints on a, Bp and (in particular they can be chosen to 
satisfy (3.3)-(3.4)), while U, and hence its asymptotic behavior, is determined by (3.7). 
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gives Q = after passing to the limit i? — > oo. The before-last equality in 
(3.10) shows that the form pBpdp+Azdz is closed, and simple-connectedness 
implies that there exists a function A such that pBpdp+ A^dz = dX, bringing 
the metric (1.1) to the form 

V + dz"^ + i.d{ip + A))^ . (3.11) 

Hence g is flat. One could now attempt to analyse (p + X near the axis of 
rotation to conclude that (p, if + X,z) provide a new global polar coordinate 
system, and deduce that g is the Euclidean metric. However, it is simpler 
to invoke the Hadamard-Cartan theorem to achieve that conclusion. 

Summarizing, we have proved: 

Theorem 3.1 Consider a metric of the form (1.1) on M = M.^, where 
{p,ip,z) are polar coordinates, with Killing vector d^, and suppose that the 
decay conditions (3.3)-(3.4) hold. If 

^R>0 

then < m < oo. Furthermore, we have m < oo if and only if 

G L^(r3) , DU, pBp^^ - A^^p G l2(m3) . 
Finally, m = if and only if g is the Euclidean metric. □ 

Remark 3.2 Theorem 2.7 shows that the coordinates required above will 
exist for a general asymptotically flat axisymmetric metric on M'^ if (2.9) 
holds with A; = 6. 

3.1 Several asymptotically flat ends 

Theorem 3.1 proves positivity of mass for axi-symmetric metrics on M^. 
More generally, one has the following: 

Theorem 3.3 Let {M,g) be a simply connected three dimensional Rieman- 
nian manifold which is the union of a compact set and of a finite number 
of asymptotic regions Mi, i = 1, . . . ,N , which are asymptotically flat in the 
sense of (2.9)-(2.10) with A; > 6. If g is invariant under an action o/U(l) 
and 

^i? > , 

then the ADM mass nii of each of the ends Mj satisfies < nii < oo, with 
mi < oo if and only if 

^R e L\Mi) , DU, pBp^, - A,^p E L^Mi) . 
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Proof: The result follows immediately from the calculations in this section 
together with Theorem 2.9: Indeed, one can integrate (3.7) on a set 

Cr := Cr \ Cyn = {-R<z<R, l/R<p<R}, 

where Cr is as in (3.6). The asymptotics (2.49) implies that the boundary 
integrals over the boundary of Ci/r give zero contribution in the limit R — > 
oo, so that (3.9) remains valid by the monotone convergence theorem in 
spite of the (mildly) singular behavior at the punctures Si of the functions 
appearing in the metric. □ 



3.2 Nondegenerate instantaneous horizons 

In order to motivate the boundary conditions in this section, recall that 
in Weyl coordinates the Schwarzschild metric takes the form (c/., e.g., [14, 
Equation (20.12)]) 



9 = 

where 



'{dp^ + dz^) 



In p — In ( 77T, sin 6 + 
1 



In 



m 



^ + + sj {z + m)2 + p2 — 2m 



\Schw 



2 



\J {z — ra)'^ + p2 + \J {^z ^ m)"^ + + 2m 
(j^Schw - m)^ - cos^ 



Schw 



2 



4y^ {z — mY + /0^\/ (-2 + + 



(3.12) 

(3.13) 
(3.14) 

(3.15) 
(3.16) 



2m + \l {^z — m)2 + + y^z + mp + p2 
In (3.13) the angle is a Schwarzschild angular variable, with the relations 



2m cos Q = (z + m)^ + — \J [z — ra)^ + p^ , 
2(rschw - m) = ^/(z + m)2 + p2 + \J {z - m)^ + p2 



/O^ = ''Schw(?^Schw - 2m) sin^ 



z = (j^Schw - m) cos 6* , 



where rgchw is the usual Schwarzschild radial variable such that e^^^^'^"' = 
1 — 2m/rschw As is well known, and in any case easily seen, C/schw is smooth 
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on 



except on the set {p = 0, —m < z < m}. From (3.13) we find, at 



fixed z in the interval —m < z < m and for small p, 



C/schw(p, z)=lnp- ln(2V(m + z)(m - z)) + 0{p^) (3.17) 

(with the error term not uniform in z). This justifies our definition: an 
interval [a, b] C £/ will be said to be a nondegenerate instantaneous horizon 
if for fixed z £ {a,b) and for small p we have 

U{p,z) = In p + U{z) +o{l) , dU{p,z) =dlnp + dU{z) + o{l) , (3.18) 

for a smooth function U. As in the Schwarzschild case the function U — a 
is assumed to be smooth across /. Thus, to compensate for the logarithmic 
singularity of U, we further assume, again for fixed z £ (a, b) and for small 
p, that there exists a function A(z) such that 



a{p, z) = U{p, z) + \{z) + o(l) . 



(3.19) 



Under those conditions the calculation of the mass formula proceeds as fol- 
lows. For k = 1, . . . , let 

h = [cfc,4] C ^ 

be pairwise disjoint intervals at which the nondegenerate instantaneous hori- 
zon boundary conditions hold. Denote by U the function, harmonic on 
\ Ufc/fc, which is the sum of Schwarzschild potentials f/schw as in (3.14), 
each with mass {d^ — Cfc)/2 and a logarithmic singularity at 1^. As in [9], 
the term \DU\'^ in (3.7) is rewritten as: 

|£»C/p = \D{U - J7 + C7)|2 = \D{U - f/)p + Di \{2U - U)D'U 

Denote by 1^ the set of points which lie a distance less than or equal to e 
from the singular set Uklk- 

Ie = {p\ d{p,UkIk) < e} • 
By inspection of the calculations so far one finds that (3.9) becomes now 



1 



-— I [DiU-U)) d'x 



H lim 

8ir e^o 



„2{a-C/)^3^ 



D'p 



D'{2U -a)- {2U - U)D'U + a— nid^^3.20) 
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In the last line of (3.20) the normal rij, taken with respect to the flat metric, 
has been chosen to point away from 7^. 

Away from the end points of the intervals Ik the logarithmic terms in U, 
U and a cancel out, leaving a contribution 



41+ / iX + f3)dz 



where \Ik\ is the length of Ik, and where we have denoted by /? the limit at 
Ufc4 oiU-U, 

kz)--= lim ^ (U{p,z)-U{p,z) 

As already pointed out, the error term in (3.17) is not uniform in z, and 
therefore it is not clear whether or not there will be a separate contribution 
from the end points of 4 to the limit as e tends to zero of the integral 
over 94. Assuming that no such contribution arises^, we conclude that the 
following formula for the mass holds: 



m 



1 

16^ 
1 

+ 8^ 



1 



D{U-U)j (fx 



(3.21) 



In the Schwarzschild case the volume integrals vanish, = 0, for z £ 
(— m, m) the function A equals 



A(.) 



1 



■In 



(m — z){z + m) 
(2m)2 



and one can check (3.21) by a direct calculation of the integral over Ii. 



3.3 Conical singularities 

So far we have assumed that the metric is smooth across the rotation axis £/. 
However, in some situations this might not be the case. One of the simplest 

*Note that this assumption, asymptotic flatness, finiteness of the volume integral in 
(3.20), and the boundary condition (3.18) on U essentially enforce the boundary condition 
(3.19) on a. 
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examples is the occurrence of conical singularities, when the regularity con- 
dition (2.45) fails to hold. It is not clear what happens with the construction 
of the coordinates (1.1) in such a case, and therefore it appears difficult to 
make general statements concerning such metrics. Nevertheless, there is at 
least one instance where conical singularities occur naturally, namely in the 
usual construction of stationary axisymmetric solutions: here one assumes 
at the outset that the space-time metric takes a form which reduces to (1.1) 
after restriction to slices of constant time; and the components of the met- 
ric are then obtained by various integrations starting from a solution of a 
harmonic map equation; c/., e.g., [8, 11, 15]. 

So consider a metric of the form (1.1) on M'^ \ {oj}, where each puncture 
di corresponds to either an asymptotically flat region or to asymptotically 
cylindrical regions (which, typically, correspond to degenerate black holes). 
Assuming that da is bounded at the axis and does not give any supplemen- 
tary contribution at the punctures, (3.9) becomes instead 



m 



1 

16^ 



+1^1 i^DUfd^x + -f 



a dz . 



„2{a-C/) ,3 



d^x 



(3.22) 



where a denotes the restriction of a to a/. 

Using (3.22) and (3.21), the reader will easily work out a mass formula 
when both conical singularities and nondegenerate instantaneous horizons 
occur. 
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